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ABSTRACT 

We have exploited the Hubble Space Telescope (HST) CANDELS J and ff-band WFC3/IR 
imaging to study the properties of (sub)millimetre galaxies within the GOODS-South field. 
After using the deep radio (VLA 1 .4 GHz) and Spitzer (IRAC 8 (jm) imaging to identify 
galaxy counterparts for the (sub)millimetre sources, we have then utilised the new CANDELS 
WFC3/IR imaging in two ways. First, the addition of new deep near-infrared photometry from 
both HST and (at AT-band) the VLT to the existing GOODS-South database has enabled us to 
derive improved photometric redshifts and stellar masses, confirming that the (sub)millimetre 
sources are massive ((-M*) = 2.2 x 10 11 ± 0.3 M Q ) galaxies at z ~ 1 — 3. Second, we have 
exploited the depth and resolution of the WFC3/IR imaging to determine the sizes and mor- 
phologies of the galaxies at rest-frame optical wavelengths A rcst > 4000A. Specifically, we 
have fitted two-dimensional axi-symmetric galaxy models to the WFC3/IR images, varying 
luminosity, axial ratio, half-light radius r 1 / 2 , and Sersic index n. Crucially, the wavelength 
and depth of the WFC3/IR imaging enables modelling of the mass-dominant galaxy, rather 
than the blue high-surface brightness features which often dominate optical (rest-frame ul- 
traviolet) images of (sub)millimetre galaxies, and can confuse visual morphological classi- 
fication. As a result of this analysis we find that > 95% of the rest-frame optical light in 
almost all of the (sub)millimetre galaxies is well-described by either a single exponential disk 
(n ~ 1), or a multiple-component system in which the dominant constituent is disk-like. 
We demonstrate that this conclusion is completely consistent with the results of recent high- 
quality ground-based A"-band imaging sampling even longer rest-frame wavelengths, and ex- 
plain why. These massive disk galaxies are reasonably extended ((^1/2) = 4.5 ± 2.5kpc; 
median r\i<x = 4.0 kpc), consistent with the sizes of other massive star-forming disks at 
z ~ 2. In many cases we find evidence of blue clumps within the sources, with the mass- 
dominant disk component becoming more significant at longer wavelengths. Finally, only a 
minority of the sources show evidence for a major galaxy-galaxy interaction. Taken together, 
these results support the view that most (sub)millimetre galaxies at z ~ 2 are simply the 
most extreme examples of normal star-forming galaxies at that era. Interestingly, the only two 
bulge-dominated galaxies are also the two lowest-redshift sources in the sample (z ~ 1), a 
result which may reflect the structural evolution of high-mass galaxies in general. 

Key words: galaxies: active - galaxies: evolution - galaxies: fundamental parameters - galax- 
ies: starburst - infrared: galaxies. 
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1 INTRODUCTION 

Since (sub)millimetre galaxies were first discovered with the 
SCUBA camera on the James Clerk Maxwell Telescope in the late 
1990s, it has been known that they are violently star-forming galax- 
ies at high redshifts. Indeed, while effective sub-mm/mm imaging 
took a long time to arrive, the strongly negative k-correction pro- 
vided by the spectral energy distribution (SED) of dust-enshrouded 
star formation (i.e. a modified black body - e.g. Hughes et al. 1993; 
Blain et al. 2002) resulted in the first sub-mm galaxies being dis- 
covered at rather high redshifts (i.e. z ~ 2 — 4; Hughes et al. 
1998). Interestingly, this meant that the supporting multi-frequency 
data available at the time was not really up to the job of revealing 
the true nature of these sources, and it has taken over a decade 
of effort with ever-improving radio surveys, optical-infrared imag- 
ing/spectroscopy with ground-based 8-10 m telescopes, and deep 
Spitzer and Hubble Space Telescope (HST) imaging to begin to clar- 
ify the physical nature of (sub)millimetre galaxies, and their role in 
the process of galaxy formation. 

This is still a work in progress. For example, it is only very 
recently that a spectroscopic redshift was finally measured for 
the brightest sub-mm galaxy discovered in the very first unbiased 
850 pm survey of the Hubble Deep Field (HDF850.1; Hughes et 
al. 1998; Downes et al. 1999; Dunlop et al. 2004; Cowie et al. 
2009) revealing it to lie at z = 5.183 (Walter et al. 2012). This 
measurement, based on CO lines, has been made possible by re- 
cent advances in sub-mm/mm spectroscopy (see also, for exam- 
ple, Riechers et al. 2010, 2012), a technique which should soon 
enable complete spectroscopic redshift distributions to be estab- 
lished for the (now substantial) samples of (sub)millimetre sources 
which have been uncovered by continuum surveys with SCUBA, 
MAMBO, LABOCA and AzTEC (e.g. Coppin et al. 2006; Greve 
et al. 2008; Weiss et al. 2009; Austermann et al. 2010). With the 
advent of the Atacama Large Millimetre Array (ALMA), high- 
resolution sub-mm/mm imaging is set to become straightforward, 
and detailed sub-mm/mm spectroscopy will soon be routine. 

However, as we now enter a new era in sub-mm/mm obser- 
vational astronomy, it is important to note that major advances 
in our knowledge of (sub)millimetre galaxies have already been 
gained from studies at optical-mid-infrared wavelengths. For ex- 
ample, aided by the improved positional accuracy provided by 
deep VLA 1.4-GHz imaging, it has already proved possible to 
identify optical/near-infrared counterparts for ~ 80% of the sub- 
mm/mm galaxies in bright samples, and hence to determine the ba- 
sic form of their redshift distribution. These studies indicate that 
(sub)millimetre galaxies peak in number density around z ~ 2 
(e.g. Chapman et al. 2005; Wardlow et al. 2011; Michalowski et 
al. 2012b; Schael et al. 2012), although there is some evidence that 
the brightest (sub)millimetre galaxies lie predominantly at z > 3 
(e.g. Smolcic et al. 2012; Koprowski et al. 2012), and the full ef- 
fects of large-scale structure have yet to be properly clarified (e.g. 
Michalowski et al. 2012b; Scott et al. 2012; Shimizu, Yoshida & 
Okamoto 2012). 

Moreover, while one might argue that redshift measure- 
ments for (sub)millimetre sources may increasingly be made at 
sub-mm/mm wavelengths, there are some crucial properties of 
(sub)millimetre-selected galaxies which can only be determined 
from deep optical-mid-infrared data. These include their stellar 
masses, formation history, and spatial distribution of their exist- 
ing stellar populations. Measurement of these physical quantities 
is of crucial importance if we are to understand the nature of 
(sub)millimetre galaxies as a function of luminosity and redshift, 



and place them in the context of the general galaxy population 
at comparable redshift/stellar-mass. Importantly, this is now an 
achievable goal, as the near/mid-infrared study of more moder- 
ate star-forming galaxies has now reached out to comparable red- 
shifts (e.g. Daddi et al. 2007; Elbaz et al. 2011) and the construc- 
tion of significant and complete mass-selected samples down to 
Mi, ~ 5 x 10 10 M Q is now possible out to z ~ 3 (e.g. Bruce 
et al. 2012). 

The most important recent advance in this context is the pro- 
vision of deep, high-quality, near-infrared imaging provided by 
WFC3/IR on HST. In particular, the CANDELS survey (Grogin 
et al. 2011; Koekemoer et al. 2011) is providing homogeneous 
WFC3/IR J125 and Hieo imaging of five fields which overlap sig- 
nificantly with existing sub-mm/mm surveys. Such imaging is of 
general importance for the study of galaxies at z ~ 1 — 3, as 
it offers high quality rest-frame optical imaging at these redshifts 
where, until recently, our knowledge of galaxy morphologies and 
colours was largely based on the rest-frame ultraviolet. However, 
it is arguably of special importance for (sub)millimetre galaxies, 
where much of the rest-frame UV light may be obscured by dust, 
and hence morphological k-corrections can be potentially dramatic. 

In this paper we exploit the CANDELS data in the GOODS- 
South field to attempt to improve our knowledge of the physi- 
cal properties of (sub)millimetre galaxies. Of the five CANDELS 
fields, GOODS-South currently offers the best overlap between 
WFC3/IR imaging and existing samples of sub-mm/mm-selected 
sources (although this situation will soon change as a result of the 
SCUBA2 Cosmology Redshift SurvejQ), 

with the added advan- 
tages of very deep Spitzer, HST ACS, and VLT imaging, and high- 
density spectroscopy (e.g. Vanzella et al. 2008). Our aim was to 
utilise the WFC3/IR imaging (Grogin et al. 2011) to first improve 
the measurement of the photometric redshifts and stellar masses 
of the known LABOCA/AzTEC sources in this field (Weiss et al. 
2009; Scott et al. 2010; Wardlow et al. 201 1; Yun et al. 2012), and 
then to determine their rest-frame optical morphologies. Crucially, 
because of the complete, homogeneous, and unbiased nature of the 
CANDELS imaging, we have also been able to compare the rest- 
frame optical morphologies of the (sub)millimetre galaxies to those 
displayed by the general population of galaxies at comparable red- 
shift and stellar mass (Bruce et al. 2012). 

This latter point is one of the key advantages offered by CAN- 
DELS, alongside, of course, the major improvement in depth and 
angular-resolution offered by WFC3/IR over previous near-infrared 
imagers. To date, deep, near-infrared imaging has only been ob- 
tained for a small number of (sub)millimetre galaxies (e.g. Smail 
et al. 2003), and generally with inadequate resolution for the reli- 
able extraction of basic morphological parameters. In an attempt 
to exploit the resolution advantage offered by HST, Swinbank et 
al. (2010) reported the results of HST NICMOS+ACS imaging of 
a sub-sample of (sub)millimetre galaxies. They did not report re- 
liable measurement of Sersic indices for individual galaxies, but 
derived median values of rn = 1.8 ± 1.0 and tij = 1.4 ± 0.8, 
and also made a measurement of Sersic index based on a stack 
of the NICMOS imaging, finding hh = 2.0 ± 0.5. Combining 
these results with measurements of half-light radii and asymme- 
try, Swinbank et al. (2010) (perhaps surprisingly, given the rather 
low median Sersic index) suggested that the stellar structure of 
(sub)millimetre galaxies is best described by a spheroid/elliptical 
galaxy light distribution. More recently, Targett et al. (2011) anal- 
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ysed deep, high-quality Jf-band images of complete subsamples of 
(sub)millimetre galaxies at redshifts z ~ 2 — 3. While the use of 
ff-band ensures imaging longward of the 4000A/Balmer break out 
to 2 ~ 4, the ground-based images have poorer resolution than 
HST data. Nevertheless, Targett et al. (2011) found the galaxies 
they studied to be reasonably extended ((r 1 / 2 ) = 3.4 ± 0.3 kpc; 
median r x /2 = 3.1 kpc) exponential disks ((n) = 1.44±0.16; me- 
dian n — 1.08). Thus, while there is improving agreement over the 
basic sizes of (sub)millimetre galaxies, there remains considerable 
(albeit possibly exaggerated) controversy over their morphologies. 

There is also an ongoing debate over the stellar masses of 
(sub)millimetre galaxies, and whether or not their extreme star- 
formation rates are triggered predominately by major galaxy- 
galaxy mergers. As discussed in Dunlop (201 1) and Michalowski et 
al. (2012a), this strikes to the heart of the nature of (sub)millimetre 
galaxies, with important implications for refining current mod- 
els of galaxy formation. For a long time it was believed that 
(sub)millimetre galaxies were the consequence of major galaxy 
mergers, in part because the brightest ULIRGs uncovered by IRAS 
in the local Universe appeared to be merger driven (e.g. Sanders & 
Mirabel 1996). However, over recent years it has become clear that 
the specific star-formation rate (s5Fi?=star-formation rate/stellar 
mass) of the general star-forming galaxy population is a factor 
~ 20 higher at z ~ 2 than at the present day (Daddi et al. 2007; 
Elbaz et al. 201 1; Karim et al. 201 1), with sSFR ~ 2 - 4 Gyr" 1 
(Gonzalez et al. 2010). Consequently, provided they have stellar 
masses M* > 10 11 M©, it can be argued that the extreme star- 
formation rates inferred for (sub)millimetre galaxies are not "unex- 
pected" at z ~ 2, and that the situation has been confused by their 
relatively early discovery. This realization has focussed increased 
attention on the accurate determination of the stellar masses of 
(sub)millimetre galaxies (e.g. Dye et al. 2008; Hainline et al. 201 1; 
Michalowski et al. 2012a) and the evidence for galaxy interactions 
(e.g. Tacconi et al. 2008; Engel et al. 2010; Alaghband-Zadeh et al. 
2012). 

Importantly, some theoretical models appear to require galaxy 
interactions, and even the additional boost provided by a top-heavy 
stellar initial mass function (IMF) to produce sufficient numbers of 
high-redshift (sub)millimetre sources (Baugh et al. 2005), but this 
may simply reflect an underlying inability of some semi-analytic 
models to produce adequate numbers of high-mass galaxies at 
z ~ 2 - 3 (Swinbank et al. 2008). Models involving AGN feed- 
back have improved in this respect (e.g. Bower et al. 2006; Cro- 
ton et al. 2006), and indeed some authors have argued that the 
(sub)millimetre galaxy population can be reconciled with the high- 
mass end of the "main sequence" of galaxy formation at z ~ 2 — 3, 
and are fed by smooth infall rather than major mergers (e.g. Dave et 
al. 2010, and see also Finlator et al. 2006; Fardal et al. 2007; Dekel 
et al. 2009). Moreover, it has been argued that major mergers are 
not sufficiently common in the redshift range of interest to explain 
the (sub)millimetre population, and that in any case the ability of 
galaxy mergers to significantly enhance the already high specific 
star-formation rates of galaxies at z ~ 2 — 3 is relatively modest, 
~ 10 — 25% (e.g. Cen 2012). Of course the true situation may 
be complex, and recently arguments have been advanced that the 
(sub)millimetre galaxy population may contain a mix of both mas- 
sive "normal" star-forming galaxies and a subset of extreme objects 
boosted by interactions (e.g. Hayward et al. 2012). 

Clearly, therefore, further observational clarification of the 
masses and morphologies of (sub)millimetre galaxies is of value, 
especially if the results can be robustly compared with the prop- 
erties of the general galaxy population. This is facilitated here by 



the work of Bruce et al. (2012), who utilised comparable CAN- 
DELS imaging to study the rest-frame optical morphologies of 
~ 200 massive (M* > 10 M galaxies in the the redshift range 
1 < z < 3. This work has shown that the CANDELS "Wide" data 
are of sufficient quality and depth to enable the robust measurement 
of the morphologies of such massive galaxies out to z ~ 3, and has 
revealed the emergence of a predominantly disk-dominated popu- 
lation by z > 2. This pattern is confirmed at somewhat lower stellar 
masses by Law et al. (2012a) who, from an HST study of 306 star- 
forming galaxies at 1.5 < z < 3.6 with M* = 10 9 - 10 11 M , 
found that star-forming galaxies at these redshifts are best de- 
scribed by an exponential (disky) surface brightness profile. How- 
ever, while there appears to be a growing consensus that massive 
galaxies at z > 2 are disk-like, there is an interesting debate over 
how these high-redshift disks compare to their local counterparts, 
since they appear to display a peaked axial-ratio distribution (e.g. 
Law et al. 2012a; Bruce et al. 2012), and have substantially higher 
velocity dispersions than low-redshift disk galaxies (e.g. Forster- 
Schreiber et al. 2009; Law et al. 2012b). There is also considerable 
debate over the prevalence and potential importance (or otherwise) 
of clumps in these high-redshift star-forming disks (e.g. Wuyts et 
al. 2012), and the extent to which they are predicted by current 
models of galaxy formation (e.g. Ceverino et al. 2012; Genel et al. 
2012). Whatever their origin and importance, it is clear that such 
high-surface brightness, and generally blue (e.g. Guo et al. 2012) 
clumps mitigate against the effective morphological classification 
of high-redshift star-forming galaxies from optical imaging (which 
samples the rest-frame near- ultraviolet at z > 1.5). 

It is clearly important to establish where (sub)millimetre 
galaxies fit into this picture. For example, as star-forming galax- 
ies they might be expected to be extended disks, and yet, based on 
their high stellar masses it is hard to avoid the conclusion that they 
are the progenitors of present-day giant elliptical galaxies. Here 
we attempt to move this subject forward using the new WFC3/IR 
CANDELS imaging to study the properties of the AzTEC and 
LABOCA-selected galaxies in the GOODS-South field. For the 
first time, we can study the rest-frame optical morphologies of 
(sub)millimetre galaxies in detail, due to the combination of high 
angular-resolution (0.1 arcsec) and depth (Hieo = 26.5; AB mag, 
5-a) provided by the CANDELS WFC3/IR imaging. In addi- 
tion, the existing rich multi-frequency data in this well-established 
field provides a significant number of spectroscopic redshifts, 
and enables the most accurate photometric redshifts achieved for 
(sub)millimetre galaxies to date (with consequently reduced uncer- 
tainties in derived stellar masses). We note that it is important to 
distinguish between the evidence for or against merger-driven star- 
bursts and the nature of the dominant underlying galaxy; in this 
paper our focus is on determining the latter and comparing the 
stellar masses and morphologies of (sub)millimetre galaxies to the 
properties of the general galaxy population at z ~ 2 — 3. We do 
also consider the imaging evidence for major galaxy-galaxy inter- 
actions, but such work is arguably better pursued with the veloc- 
ity information gleaned from optical and/or mm spectroscopy. It 
does, however, need to be pursued in an unbiased manner includ- 
ing both (sub)millimetre bright and non (sub)millimetre-detected 
massive galaxies at z ~ 2, before definitive conclusions can be 
drawn regarding the importance of major mergers in producing the 
(sub)millimetre galaxy population. 

This paper is structured as follows. In Section 2 we summarise 
the available (sub)millimetre survey data within GOODS-South, 
the new CANDELS HST imaging, and the additional key support- 
ing data available within this field. Next, in Section 3, we explain 
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how we determined the the radio/IRAC and hence WFC3/IR H- 
band galaxy counterparts for the (sub)millimetre sources. Then, in 
Section 4 we describe the two-dimensional modelling used to ex- 
tract the basic morphological properties of the galaxies, and the 
SED fitting to the multi-frequency data used to determine redshifts 
and stellar masses. We present the results of our analysis and, in 
Section 5, also attempt to place our findings in the context of studies 
of other galaxy populations at both high and low redshift. Finally, 
our main conclusions are summarised in Section 6. Throughout we 
quote magnitudes in the AB system (Oke 1974), and assume a cos- 
mological model with Hq = 70 kms^Mpc -1 , S7a = 0.7 and 

n m = 0.3. 



2 DATA 

2.1 (Sub)millimetre data 

The full 30 x 30 arcmin ECDFS has been mapped at a wavelength of 
870 by the LABOCA ECDFS (sub)millimetre Survey (LESS) 
(Weiss et al. 2009). The 12-m diameter of the APEX telescope 
(Gusten et al. 2006) delivers a 19.2 arcsec FWHM beam at this 
wavelength. The LESS image has a uniform depth of as70 — 

1.2 mjy beam -1 (as measured in the map, and hence including 
confusion noise). 

A 512 arcmin 2 1.1 mm map of GOODS-South has now also 
been obtained with the AzTEC camera on the ASTE, which deliv- 
ers a 30 arcsec FWHM beam at this wavelength (Scott et al. 2010). 
Within the deepest central 270 arcmin region of the AzTEC map, 
the rms noise ranges from 0.48 to 0.73 mjy beam -1 , making this 
GOODS-South millimetre map the deepest contiguous region ever 
mapped at 1 . 1 mm. 

The (sub)millimetre galaxy sample considered here consti- 
tutes the subset of published 870 /im (LABOCA) and 1.1mm 
(AzTEC) sources which lie within the areas covered by the 
HST WFC3/IR imaging of the CANDELS deep (10-orbit depth), 
wide (2-orbit depth), and ERS (2-orbit depth) surveys (total area 
~ 164 arcmin 2 ) centred on the southern field of the Great Obser- 
vatories Origins Deep Survey (GOODS; Dickinson et al. 2004) at 
RA03 h 32 m 30 s , Dec - 28° 48' (J2000). 

Our final GOODS-South (sub)millimetre sample comprises 
the 24 sources listed in Table 1 (AzTEC. GS 19 is tabulated twice 
here because it has two alternative IDs, as explained in Section 
3). Because the AzTEC imaging is slightly deeper, only 3 of the 
sources are uniquely supplied by LESS (as indicated in Table 1). 
There are, however, 6 sources in common, and here the LESS imag- 
ing is still of value as it supplies a somewhat better position due to 
the smaller beam size of the shorter-wavelength survey. 



2.2 CANDELS data 

The Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey 
(CANDELS; Grogin et al. 201 1; Koekemoer et al. 201 1) is an HST 
Multi-Cycle Treasury Program designed to image portions of five 
legacy fields (GOODS-N, GOODS-S, COSMOS, UDS, and EGS) 
with the Wide Field Camera 3 (WFC3) in the near infrared. The 
CANDELS data have been reduced and drizzled to a 0.06 arcsec 
pixel scale (Koekemoer et al. 2011) to create the high-resolution 
mosaics. Despite the variation in exposure times across the field, 
even the shallowest region reaches a 5-a AB limiting magnitude 
of 26.5 in the 7Ji6o-band, two magnitudes deeper than the existing 



ground-based ISAAC coverage at the same wavelength (Retzlaff et 
al. 2010). 

2.3 Supporting multi-frequency data 

2.3.1 Radio: VLA 1.4-GHz imaging 

As demonstrated by the follow-up of the sub-mm/mm surveys 
undertaken with SCUBA, IRAM, LABOCA and AzTEC, very 
deep VLA imaging is necessary (and frequently sufficient) for 
the successful identification of the galaxy counterparts of sources 
detected with the large beams delivered by single-dish far- 
infrared/(sub)millimetre facilities (e.g. Ivison et al. 2002, 2007; 
Michalowski et al. 2012b; Yun et al. 2012). This works for three 
reasons. First, star-forming galaxies produce copious quantities of 
synchrotron emission. Second, even in the deepest available radio 
maps, 1.4 GHz sources still have sufficiently low surface-density 
that positional coincidences within a 'reasonable' search radius are 
generally statistically rare. Third, if a robust radio counterpart is 
found, the ~ 1 arcsec positional accuracy provided by the VLA at 

1.4 GHz essentially always yields an unambiguous optical/infrared 
galaxy counterpart for further study. 

Very deep (cri.4 ^ 7.5 /iJy), high-resolution 1.4 GHz imag- 
ing of GOODS-South is available at the centre of the ECDFS as 
described by Kellermann et al. (2008) and Miller et al. (2008). We 
adopt the radio catalogue from Dunlop et al. (2010) who, using 
the techniques described by Ibar et al. (2009), created a radio cata- 
logue down to a 4-a limit of 30 fi Jy (at which depth the cumulative 
source density on the sky is ~ 0.8 arcmin -2 ). This catalogue was 
then searched for potential radio counterparts to the sub-mm/mm 
sources using the method described in Section 3. We note that even 
deeper VLA 1.4 GHz imaging in GOODS-South has recently been 
presented by Yun et al. (2012), who have also searched for the 
counterparts to the AzTEC sources discussed here (using the same 
statistical technique). Where appropriate, we compare and cross- 
check our galaxy identifications in Section 3. 

2.3.2 Optical: HST ACS imaging 

Deep optical imaging over the GOODS-South field was ob- 
tained (prior to CANDELS) with the Advanced Camera for Sur- 
veys (ACS) on board HST in 4 different filters: F435W (S435), 
F606W (Veoe), F775W (i 775 ) and F850LP (z 850 ). These data 
were taken as part of GOODS (Giavalisco et al. 2004) and reach 
5-a limiting (AB) magnitudes (within a 1 arcsec diameter aperture) 
of 26.9, 26.9, 26.1, and 25.8, respectively. 

2.3.3 Ground-based K s -band VLT imaging 

Ground-based near-infrared imaging of almost all of the CAN- 
DELS GOODS-South field was completed several years ago with 
the ISAAC camera on the VLT (Retzlaff et al. 2010). These 
data cover ~ 143 arcmin 2 within the ~ 164 arcmin 2 CANDELS 
GOODS-South field. Although ground-based, these data remain of 
value because the HST cannot image longward of the if -band. In 
addition, the image quality of the ISAAC data is excellent, with a 
median FWHM PSF value of ~ 0.5 arcsec (variable from ~ 0.37 
to ~ 0.7 arcsec across the field; Bouwens et al. 2008). Partly as a 
result of this good image quality, the point-source sensitivity of this 
imaging is deep, reaching a 5-a detection level of K s ~ 24.7 (AB 
mag). 

As part of the larger CANDELS project, even deeper 
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Table 1. Positional information on the 24-source combined AzTEC+LABOCA GOODS-South (sub)millimetre galaxy sample utilised in this paper, and the 
adopted galaxy counterparts in the HST CANDELS Hiqo imaging. The first column gives the name of the AzTEC source from Scott et al. (2010), with the 
next two columns giving the AzTEC position (AzTEC. GS 19 is tabulated twice here because it has two alternative IDs, as explained in Section 3). Column 4 
then gives the source name from the LESS survey (Weiss et al. 2009), where here the name itself contains the original LABOCA survey position. Columns 5 
and 6 then give the positions of the Hu;o galaxy counterpart we have selected for morphological study. These galaxy identifications were selected as described 
in Section 3. Column 7 then gives the probability that the adopted galaxy counterpart is merely a chance association, while column 8 indicates whether the 
galaxy identification was selected on the basis of the radio or Spitzer imaging, and the source of the quoted P value; (a) our identification process, or (b) from 
Yunetal. (2012). 



AzTEC 


RA 


Dec. 


LABOCA 


CANDELS 


CANDELS 


P 


Type 


ID 


(J2UUU) 


(J20U0) 


ID 


KA (J20U0) 


Dec (J 2000) 








AzTEC.GS03 


03:32:47.86 


-27:54:19.3 


LESSJ033248. 1-275414 


03:32:47.992 


-27:54:16.42 


0.017 


a - 


1.4 GHz 


AzTEC.GS06 


03:32:25.73 


-27:52:19.4 


LESSJ033225.7-275228 


03:32:25.258 


-27:52:30.51 


0.062 


a - 


1.4 GHz 


AzTEC.GS08 


03:32:05.12 


-27:46:45.8 


LESSJ033205. 1-274652 


03:32:04.873 


-27:46:47.37 


0.015 


a - 


1.4 GHz 


AzTEC.GSll 


03:32:15.79 


-27:50:36.8 




03:32:15.319 


-27:50:37.12 


0.067 


a - 


1.4 GHz 


AzTEC.GS12 


03:32:29.13 


-27:56:13.8 


LESSJ033229.3-275619 


03:32:29.290 


-27:56:19.47 


0.032 


b - 


1.4 GHz 


AzTEC.GS13 


03:32:11.91 


-27:46:16.9 




03:32:11.908 


-27:46:15.51 


0.009 


b - 


1.4 GHz 


AzTEC.GS16 


03:32:37.67 


-27:44:01.8 




03:32:38.006 


-27:44:00.57 


0.066 


a - 


1.4 GHz 


AzTEC.GS17 


03:32:22.31 


-27:48:16.4 




03:32:22.566 


-27:48:14.83 


0.026 


a/b 


- 8/24 /an 


AzTEC.GS18 


03:32:43.58 


-27:46:36.9 


LESSJ033243.6-274644 


03:32:43.536 


-27:46:38.90 


0.035 


b- 


1.4 GHz 


AzTEC.GS19-l 


03:32:23.21 


-27:41:28.8 




03:32:22.879 


-27:41:24.96 


0.074 


a - 


1.4 GHz 


AzTEC.GS19-2 


03:32:23.21 


-27:41:28.8 




03:32:22.708 


-27:41:26.39 


0.088 


a - 


1.4 GHz 


AzTEC.GS21 


03:32:47.60 


-27:44:49.3 




03:32:47.592 


-27:44:52.26 


0.026 


a - 


1.4 GHz 


AzTEC.GS22 


03:32:12.60 


-27:42:57.9 




03:32:12.545 


-27:43:05.99 


0.116 


b- 


1.4 GHz 


AzTEC.GS23 


03:32:21.37 


-27:56:28.1 


LESSJ033221.3-275623 


03:32:21.574 


-27:56:23.94 


0.054 


a - 


1.4 GHz 


AzTEC.GS24 


03:32:34.76 


-27:49:43.1 




03:32:34.163 


-27:49:39.52 


0.141 


b- 


1.4 GHz 


AzTEC.GS26 


03:32:15.79 


-27:43:36.6 




03:32:16.281 


-27:43:43.41 


0.624 


a - 


8 (im 


AzTEC.GS27 


03:32:42.42 


-27:41:51.9 




03:32:41.632 


-27:41:51.45 


0.623 


a - 


8 fim 


AzTEC.GS28 


03:32:42.71 


-27:52:06.8 




03:32:41.853 


-27:52:02.47 


0.734 


a - 


8 (im 


AzTEC.GS30 


03:32:20.94 


-27:42:40.8 




03:32:20.658 


-27:42:34.40 


0.082 


b- 


1.4 GHz 


AzTEC.GS34 


03:32:29.77 


-27:43:13.1 




03:32:29.479 


-27:43:22.10 


0.171 


a - 


1.4 GHz 


AzTEC.GS35 


03:32:26.90 


-27:40:52.1 




03:32:27.181 


-27:40:51.42 


0.008 


b- 


1.4 GHz 


AzTEC.GS38 


03:32:09.26 


-27:42:45.5 




03:32:09.705 


-27:42:48.10 


0.015 


a - 


1.4 GHz 








LESSJ033217.6-275230 


03:32:17.617 


-27:52:28.52 


0.016 


a - 


8 /an 








LESSJ033219.0-275219 


03:32:19.032 


-27:52:13.97 


0.001 


b- 


1.4 GHz 








LESSJ033243.3-275517 


03:32:43.179 


-27:55:14.49 


0.019 


a - 


1.4 GHz 



ground-based if s -band imaging of GOODS-South (covering ~ 
88arcmin 2 , within the CANDELS imaging region) has recently 
been obtained with the HAWK-I imager on the VLT. These data 
have been reduced and processed as described in Fontana et al. 
(2012). Again taken in excellent seeing (~ 0.45 arcsec), this imag- 
ing reaches a 5-a limiting point-source depth of K s ~ 26.0 (AB 
mag), comparable to the CANDELS HST /fieo-band data. 



2.3.4 Mid-infrared: Spitzer imaging 

Again as part of GOODS, ultra-deep Spitzer imaging with the In- 
frared Array Camera (IRAC: Fazio et al. 2004) has been obtained 
over the whole of the CANDELS GOODS-South field, in all 4 
IRAC channels (3.6, 4.5, 5.6, and 8.0 ^m). The IRAC 5-a detection 
limits are S 3 .6 ^ 25.9, 5*4.5 ^ 25.5, S 5 . 6 ~ 23.3, S 8 .o ^ 22.9 
(AB mag). 

The 24 fim Spitzer MIPS data originally obtained as part of 
GOODS has been augmented and incorporated within the Spitzer 
Far-Infrared Deep Extragalactic Legacy (FIDEL0 survey (Mag- 
nelli et al. 2009), and reach 5-a detection limits of S24 — 30 /iJy. 



2 PI M. Dickinson, see |http : / /www . noao ."ed u/noao/ f idel/ 



3 SELECTION OF GALAXY COUNTERPARTS 
3.1 Radio identifications 

For the reasons given in subsection 2.3.1, the primary method 
adopted to achieve the improved positional information required 
to connect the (sub)millimetre sources to unique galaxy counter- 
parts in the optical/near-infrared imaging is via the identification 
of statistically-significant counterparts in the deep radio 1.4 GHz 
imaging. Here this is even more challenging than in the surveys 
undertaken with SCUBA on the JCMT (e.g. Ivison et al. 2007) 
because of the larger size of the AzTEC and LABOCA beams 
in the GOODS-South surveys (FWHM=30 arcsec at 1.1mm and 
FWHM=19.2 arcsec at 870 respectively). 

However, the radio data described above in subsection 2.3.1 
are certainly deep enough to find counterparts to the vast majority 
of the (sub)millimetre sources, especially when supplemented by 
the even deeper radio imaging of Yun et al. (2012), and the only 
real issue is establishing an acceptable level of statistical security. 

The anticipated uncertainty in the (sub)millimetre positions is 
given by 

a po3 {= Aq = AS) = — 
where 9 is the FWHM of the (sub)millimetre beam, and S/N is 
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AzTEC.GS08 




Figure 1. Examples of the results of two-dimensional modelling of what appears to be a single-component clumpy disk galaxy (AzTEC.GS08; upper row) 
and a multi-component system (LESSJ033243; lower row). The left-hand panels show colour images (6" x6") based on the I + V , J125 and -H160 images. 
The next column shows grey-scale Higo-band CANDELS postage-stamp images (12" X 12") centred on the (sub)millimetre galaxy counterparts. The third 
column shows the best-fitting two-dimensional models. The right-hand panels show the residual images after subtraction of the best-fitting models from the 
data. The panels in columns 2, 3 and 4 are shown on the same greyscale. 



the statistically-deboosted signal:noise ratio of the (sub)millimetre 
source. To be conservative, we adopted aS/N value of 3, and hence 
searched for counterparts with a 2.5-cr search radius of 15arcsec 
around the AzTEC sources and lOarcsec around the LABOCA 
sources (a radius expected to contain 95% of all genuine radio 
counterparts). 

Having located all potential radio counterparts within the 
search radii, we then calculated the probability that each poten- 
tial counterpart could have been found simply by chance using the 
method described by Downes et al. (1986) (see also Dunlop et al. 
1989). This involves first calculating the raw Poisson probability 
that a radio source of the observed 1.4 GHz flux density would 
be discovered by chance at the measured distance from the nom- 
inal (sub)millimetre source position, and then correcting this prob- 
ability for the number of ways such a statistical coincidence could 
have been discovered given the available search parameter space 
defined by the maximum search radii, the radio-source number 
density at the limiting search flux-density (N(Sia > 30/iJy) ~ 
0.8 arcrnin -2 ), and the form of the radio source counts over the 
flux-density range of interest (here we adopt a power-law index 
of 1.4). For the deep data under study here, this correction is of- 
ten substantial, typically increasing P by a factor of a few. This 
technique has been applied previously to estimate the robustness 
of radio identifications for SCUBA sources (e.g. Ivison et al. 2002, 
2007), and BLAST sources (Dunlop et al. 2010), and is also the 
technique adopted by Yun et al. (2012) in their own work on the 
radio identification of the GOODS-South AzTEC sources. 

The calculated value of P is the probability that the observed 
association is the result of chance. Even a very low value of P does 
not prove that the radio source is the (sub)millimetre source. Never- 
theless, a low value of P clearly does imply that the radio source is 
likely related to the(sub)millimetre source in some way. This could 



be true for several reasons. First, the radio source could be the sin- 
gle true counterpart of the (sub)millimetre source. Alternatively it 
could be one of a group of two or more galaxies which contribute to 
the (sub)millimetre flux-density peak (of particular relevance here, 
given the large beam-sizes). Third, the statistical result could be a 
consequence of some secondary association (e.g. clustering with 
the true (sub)millimetre sources, or the result of gravitational lens- 
ing). Despite these potential concerns over statistical interpretation, 
we re-emphasise that our search for radio counterparts is not mo- 
tivated purely by their statistical rarity and good positional accu- 
racy, but also by the physical evidence that dust-enshrouded star- 
forming galaxies also produce copious quantities of synchrotron 
emission (resulting in the well-known far-infrared:radio luminosity 
correlation - e.g. Ivison et al. 2010). A sensible hypothesis, there- 
fore, is that an apparently statistically secure radio counterpart to 
a (sub)millimetre source is at least a significant contributor to the 
observed (sub)millimetre emission. 

As summarised in the final two columns of Table 1, the radio 
imaging yields counterparts for 19 out of the 24 (sub)millimetre 
sources, all but 3 of which have P < 0.1. The sum of all the P 
values for the 19 radio identifications suggests, statistically, that 
only one of these sources may have been mis-identified. 



3.2 Spitzer 8/im identifications 

This leaves five sources which lack a galaxy identification. To 
identify the most probable counterparts for these we follow re- 
cent work by Michalowski et al. (2012b) and Koprowski et al. 
(2012) which confirms that the brightest Spitzer 8pm source within 
the search radii (as defined above) can be reliably adopted as the 
correct galaxy counterpart, even when a radio source (or indeed 
a 24 pm counterpart) is absent (possibly because the source lies 
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at a high redshift where the k-correction from the steep-spectrum 
synchrotron spectrum has rendered the radio counterpart unde- 
tectable). This appears to work because sub-mm/mm galaxies are 
among the most massive at any epoch z > 1 and even at, say, 
z ~ 4, the observed 8 fim band continues to sample the rest-frame 
SED of a galaxy close to the peak emission of its mass-dominant 
stellar population (i.e. at \ res t — 1.6 /iin). 

Indeed for 18 out of the 19 sources in Table 1 which have ra- 
dio identifications as described above, selecting the brightest 8 /im 
counterpart would have delivered the same optical/near-infrared 
galaxy identification, and the only reason we have listed these ob- 
jects as primarily radio identifications is that the P values are lower 
(because the surface density of radio sources at the limit of the ra- 
dio imaging is significantly lower than that of the 8 pm sources at 
the limit of the (deep) GOODS-South Spitzer imaging). 

We thus searched for, and successfully located, 8 (im counter- 
parts for the remaining five sources. Two of these are statistically 
secure, with P < 0.05, but the remaining three have rather high P 
values, and must be regarded as only tentative identifications. 

3.3 Final galaxy sample 

Table 1 thus summarises the 25 proposed identifications to the 24 
(sub)millimetre sources which are the subject of the further mor- 
phological and SED analyses presented in this paper. The reason 
these numbers do not match is simply that, as indicated in Table 1, 
we have chosen to retain the two alternative galaxy counterparts for 
AzTEC.GS19 because they are both equally probable, and at least 
one such source blend is not unexpected in this sample given the 
large size of the AzTEC beam. Of our 24 (sub)millimetre sources, 
18/24 have secure galaxy identifications (P < 0.1) while 6/24 
are statistically insecure (P > 0.1). While it would be desirable 
to insist on a significance threshold of P < 0.05, we here adopt 
P < 0.1 as appropriate given the large search radius required by 
the AzTEC and LABOCA beams (12 sources still have P < 0.05). 
Despite the fact that a sum of the P values indicates that at most 
only 3 sources may have been misidentified, the less secure iden- 
tifications (P > 0.1) are explicitly flagged as a subset in all the 
subsequent analysis presented in this paper (it transpires that our 
results are unaffected by their inclusion or exclusion). 

Of the galaxy identifications, two objects, AzTEC. GS 12 and 
AzTEC. GS28, were found to be unresolved in the //-band imag- 
ing, rendering them unsuitable for two-dimensional modelling, 
although we retained them in the sample for SED fitting (note 
that AzTEC. GS28 has the formally least secure galaxy identifica- 
tion in our whole sample, so may well have been misidentified; 
AzTEC.GS12 is presumably an AGN, as it has a secure radio coun- 
terpart). Finally we note that AzTEC. GS03, although a secure ra- 
dio identification, appears too faint in the //-band (and all shorter 
wavelengths) for morphological modelling, although its solid IRAC 
detections still allowed meaningful SED fitting. 

We are thus left with a sample of 21 (sub)millimetre sources, 
with 22 galaxy counterparts suitable for two-dimensional mod- 
elling as described in the next section. 



4 MODEL FITTING 

4.1 Two-dimensional modelling fitting 

Morphological parameters for the (sub)millimetre galaxy sample 
were obtained from the GOODS-South deep CANDELS PT-band 



data using GALFIT (Peng et al. 2002, 2010). The GALFIT code 
employs \ 2 minimisation to create a best-fit numerical model, 
where free parameters are tuned to match the light distribution of 
a source. For comparison with the data the model obviously needs 
to be convolved with the appropriate imaging point spread function 
(PSF); here the PSF for the CANDELS imaging was derived using 
the IRAF package PSF, fitting several stars from the science image 
with multiple statistical distribution types to determine the best fit- 
ting model. The resulting models were then centroided, averaged, 
and processed with look-up tables to reproduce any asymmetries. 

The results of the two-dimensional modelling for all 22 vi- 
able (sub)millimetre galaxy counterparts are presented in Table 2, 
with examples of models and model-subtracted residuals for sin- 
gle and multi-component fits shown in Fig. 1 (see Fig. Al in the 
Appendix for the model fits to all objects). We find that > 95% 
of the rest-frame optical light in almost all of the (sub)millimetre 
galaxies is best-described by either a single-component exponen- 
tial disk (10/22 objects), or multiple component systems where the 
dominant constituent is disk-like (10/22). 

For comparison with existing studies of larger samples of mas- 
sive galaxies, and as a test of Sersic index measurements obtained 
from lower-resolution data (i.e. Targett et al. 2011), all multiple- 
component sources were also fitted with a single-component model 
forced to attempt to reproduce the entire source (see SCP entries in 
Table 2). We find that both the average size and Sersic index of 
these single component fits match those of the dominant compo- 
nent in the multi-object fits to within a typical accuracy of 5%. This 
similarity highlights the dominant nature of the disk-like compo- 
nent in the multiple component systems. It also suggests that single 
Sersic fits to complex systems in moderate-resolution data (e.g. the 
~ 0.4arcsec seeing ground-based imaging analysed by Targett et 
al. 2011), while missing the fine-structure seen at HST resolution, 
can successfully recover the basic morphology of the underlying 
mass-dominant galaxy (see the further analysis of this issue pre- 
sented in Section 5.1). 

Among the 10 multiple-component systems we find an even 
division between objects whose secondary components display a 
Sersic index consistent with a disky morphology (possibly indi- 
cating a companion, possibly interacting or merging galaxy), and 
objects with secondary components which appear to be best de- 
scribed by extremely low Sersic indices outside the normal range 
displayed by galaxies. The anomalous Sersic index of these sec- 
ondary components could be indicative of regions of intense star- 
formation or giant clumps believed to arise in star-forming galaxies 
at high-redshift (Ceverino et al. 2012, Hopkins et al. 2012). This 
separation into galactic and non-galactic secondary components is 
not possible via visual classification, and demonstrates another ad- 
vantage of performing detailed multiple-component variable-Sersic 
index model fitting to complex objects. 

In order to further check the accuracy and reliability of model 
parameter recovery, we have also measured the Petrosian radius and 
ri/2 via a curve-of-growth analysis directly from the Pieo-band 
data. The Petrosian radius was defined as the point at which the ra- 
tio of the surface brightness in an annulus at radius r to the mean 
surface brightness within an aperture of radius r equals 0.2 (Blan- 
ton et al. 2001, Yasuda et al. 2001). It was found that the values of 
ri/2 as determined from the two-dimensional modelling, curve-of- 
growth, and Petrosian techniques were consistent to within 10%. 
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Table 2. Results from two-dimensional modelling of the CANDELS Higo galaxy counterparts (including sub-components where significant) to the AzTEC 
and LABOCA selected (sub)millimetre galaxies in CANDELS GOODS-South. All multiple-component sources were also fitted with a single-component 
model forced to attempt to reproduce the entire source (entries labelled SCP, shown just below the results from the separate components of the source name 
listed above). Column 1 gives the source name. Column 2 lists the total model host-integrated i/igo-band magnitude (AB). Column 3 lists the semi-major axis 
scalelength (half-light radius) of the host galaxy fit in kpc. Column 4 lists the value of the Sersic index (rt). Column 5 gives the axial ratio of the host galaxy. 
Column 6 gives the value of reduced \1 ror eacn model fit, where multiple-component fits share a common xt f° r the entire combined fit. 



Name 


H 160 / AB mag 


r 1/2 /kpc 


Sersic n 


Axial ratio 


2 


AzTEC.GS06 


23.64 


3.5 


0.8 


0.34 


0.997 




24.49 


11.8 


9.4 


0.72 


0.997 


SCP 


23.76 


3.3 


0.7 


0.76 


1.030 


AzTEC.GS08 


21.92 


6.3 


1.0 


0.51 


1.196 


AzTEC.GSll 


23.91 


3.8 


0.7 


0.48 


1.270 


AzTEC.GS12 


unresolved 










AzTEC.GSB 


24.13 


3.1 


1.2 


0.79 


0.885 




24.81 


3.2 


0.8 


0.59 


0.885 


SCP 


23.76 


4.6 


0.7 


0.53 


0.916 


AzTEC.GS16 


24.39 


4.6 


1.0 


0.93 


0.954 


AzTEC.GS17 


23.31 


1.6 


1.4 


0.96 


1.150 


AzTEC.GS18 


24.38 


5.3 


1.1 


0.49 


0.962 


AzTEC.GS19-l 


22.18 


4.1 


0.9 


0.72 


1.717 


AzTEC.GS19-2 


23.15 


2.8 


1.2 


0.65 


1.717 


AzTEC.GS21 


22.94 


2.6 


1.3 


0.65 


2.048 




24.25 


3.2 


0.3 


0.56 


2.048 




24.50 


2.3 


0.2 


0.34 


2.048 


SCP 


22.42 


3.7 


1.5 


0.82 


2.635 


AzTEC.GS22 


23.70 


4.5 


1.2 


0.40 


1.036 




24.22 


3.5 


0.3 


0.44 


1.036 


SCP 


23.32 


6.0 


0.7 


0.27 


1.065 


AzTEC.GS23 


22.67 


4.6 


1.3 


0.67 


1.193 




23.16 


1.7 


0.2 


0.79 


1.593 




23.28 


3.8 


0.7 


0.31 


1.593 




23.38 


3.0 


1.0 


0.62 


1.593 


SCP 


21.54 


5.9 


1.1 


0.67 


1.558 


AzTEC.GS24 


24.24 


1.7 


1.0 


0.59 


1.913 


AzTEC.GS26 


23.01 


6.6 


1.3 


0.30 


1.148 


AzTEC.GS27 


20.64 


3.7 


4.0 


0.94 


1.875 




21.70 


7.0 


1.0 


0.73 


1.875 


AzTEC.GS28 


unresolved 










AzTEC.GS3() 


21.93 


6.3 


1.4 


0.66 


1.457 




24.51 


3.5 


0.2 


0.22 


1.457 


SCP 


21.88 


5.4 


1.1 


0.68 


1.824 


AzTEC.GS34 


21.92 


1.6 


1.9 


0.84 


1.193 


AzTEC.GS35 


23.69 


1.8 


1.3 


0.48 


1.426 




23.76 


0.8 


0.9 


0.55 


1.426 


SCP 


23.05 


2.8 


0.7 


0.81 


2.789 


AzTEC.GS38 


19.32 


4.4 


4.0 


0.95 


2.219 




21.46 


7.2 


1.0 


0.62 


2.219 


LESSJ033217 


20.16 


5.9 


1.0 


0.82 


1.274 




22.65 


1.4 


0.3 


0.28 


1.274 


SCP 


19.88 


6.4 


2.0 


0.78 


1.489 


LESSJ033219 


24.31 


2.4 


1.4 


0.48 


1.058 




24.78 


1.1 


1.3 


0.48 


1.058 


LESSJ033243 


22.25 


7.6 


1.1 


0.53 


1.120 




23.47 


1.6 


0.8 


0.66 


1.120 




24.17 


1.6 


0.4 


0.38 


1.120 




24.37 


1.1 


0.4 


0.38 


1.120 


SCP 


22.04 


7.5 


0.7 


0.38 


2.686 



4.2 SED fitting 

The extensive and deep multi-frequency data in the GOODS-South 
field, now further enriched with the deep CANDELS near-infrared 
imaging, allows high-quality SED fits to the (sub)millimetre galaxy 



counterparts to determine redshifts, stellar masses, and other model 
parameters. The photometry for the SED fitting was extracted 
via 2 arcsec diameter aperture photometry at the position of each 
CANDELS fiieo counterpart in the HST ACS ^435^ (£435), 
F606W (Veoe), F775W (1775) and F850LP (z 8B o) bands, the 
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Figure 2. Example plots for multiple-component SED fitting to LESSJ033219 (top left), LESSJ033243 (top right), AzTEC.GS08 (bottom left), and 
AzTEC.GS13 (bottom right). In each plot, the red line represents the contribution from the older stellar population, while the blue line shows the contri- 
bution from the most recent burst of star formation. For each source the secondary panel shows the value of x 2 as a function of redshift (marginalized over all 
other free parameters). 



CANDELS F125W J-band and F160W #-band, and in all 4 
IRAC channels (3.6, 4.5, 5.6, and 8.0 /im). The 2arcsec aperture 
did not require any significant aperture correction in the optical and 
near-infrared HST imaging, while the IRAC data were de-confused 
to comparable resolution using the CANDELS i/iso image as a 
template (as described in McLure et al. 201 1). 

The SED fitting procedure we then applied to derive photo- 
metric redshifts and stellar masses is based largely on the public 
package HYPERZ (Bolzonella, Miralles, & Pello 2000). The ob- 
served photometry was fitted with synthetic galaxy templates gen- 
erated with the Bruzual & Chariot (2003) stellar population models. 
We used a variety of star-formation histories: single and double- 
component instantaneous bursts, and exponentially-declining star- 
formation. Following the analysis of (sub)millimetre galaxy coun- 
terparts presented in Michalowski et al. (2012a), we assumed solar 
metallicity and a Chabrier (2003) initial mass function (IMF). At 
each redshift step we allowed galaxy templates with ages ranging 
from 1 Myr to the age of the Universe at that redshift. For dust 



reddening, we adopted the prescription from Calzetti et al. (2000), 
within the range < Ay < 2, and included absorption due to HI 
clouds along the line-of-sight according to Madau et al. (1996). The 
results of the multi-frequency SED fitting for the (sub)millimetre 
galaxies are presented in Table 3, with examples of the SED fits 
shown in Fig. 2, including AzTEC.GS08 and LESSJ033243, the 
two example objects already presented in Fig. 1. 

In Table 3 we also give the spectroscopic redshifts for those 
sources for which robust spectroscopic results are available in the 
literature. Due the intensive nature of the spectroscopic campaigns 
in GOODS-South, almost half of our (sub)millimetre sample (i.e. 
10 sources) have spectroscopic redshifts (an unusually good situa- 
tion for a (sub)millimetre galaxy sample). In Fig. 3 we plot z p hot 
v z spec for these 10 sources, which serves to demonstrate the qual- 
ity of the photometric redshifts derived from the unusually high- 
quality multi-frequency photometry now available in the GOODS- 
South field. Clearly this figure provides confidence in the photo- 
metric redshifts derived for the remaining sources. 



© 0000 RAS, MNRAS 000, 000-000 



10 T.A. Targett, et al. 




12 3 4 5 

z / spectroscopic 



Figure 3. A comparison of the photometric and spectroscopic redshifts for 
the 10 sources in our sample which already possess good-quality optical 
spectroscopy (from the literature sources indicated in Table 3). The error 
bars represent the 1-ct uncertainties in z phot as derived from the SED fits, 
as given in Table 3. 

5 DISCUSSION 

5.1 Size and Morphology 

We find that > 95% of the rest-frame optical light in al- 
most all of the (sub)millimetre galaxies is well-described by 
either a single-component exponential disk (10/22 objects), or 
multiple-component systems where the dominant constituent is 
disk-like (10/22 objects). This leaves two sources (AzTEC.GS27 
and AzTEC.GS38) which are best fitted by a disk+bulge model 
with a dominant bulge component, but interestingly one of these 
(AzTEC.GS27) is one of the 3 least secure identifications in the 
sample, while the other (AzTEC.GS38) is the only galaxy in the 
sample with a confirmed redshift z < 1. In summary, all of the 
(sub)millimetre galaxies in our sample at z > 1.5 appear to be 
disks. 

To place these results in context, the distribution of Sersic in- 
dices for our (sub)millimetre galaxies is plotted in Fig. 4 along with 
the distribution of Sersic indices derived by Bruce et al. (2012) 
for a sample of 83 comparably-massive (M* > lO n M ) galax- 
ies at 2 ~ 2 from _Hi6o CANDELS imaging of essentially identi- 
cal quality and depth. Bruce et al. (2012) have already noted that 
a significant fraction of massive galaxies at these redshifts appear 
to be disk-dominated, and that a substantial majority of these are 
star-forming. However, this figure shows that the (sub)millimetre 
selected galaxies studied here are particularly disk-like in the sense 
that their Sersic indices lie very close to unity. Indeed, it appears 
that (sub)millimetre selection is a good way of pre-selecting a sub- 
stantial fraction of the most disk-dominated massive galaxies at 
these redshifts. 



Comparing these new results to those already in the litera- 



Table 3. The results of SED fitting to the optical-mid-infrared photometry 
of the AzTEC and LABOCA selected (sub)millimetre selected galaxies in 
CANDELS GOODS-South. Column 1 gives the source name. Columns 2, 
3, 4, and 5 list the photometric redshifts, the l-cr uncertainty in the pho- 
tometric redshifts, the stellar masses, and \ 2 f° r the best-fitting SED tem- 
plate. Column 6 gives the spectroscopic redshift where available, and the 
reference from which this redshift was obtained: a from Vanzella et al. 
(2008), b from Huang et al. (2012), c from Balestra et al. (2010), d from 
Straughn et al. (2009), e from Kurk et al. (2009), and f from Kriek et al. 
(2007). Since spectroscopic redshifts are available for 10 sources (i.e. es- 
sentially half the sample), we show a plot of z spec versus 2 p hot f° r these 
objects in Fig. 3, which demonstrates the accuracy of the photometric red- 
shifts as derived from the available high-quality photometry. 
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AzTEC.GS08 


2.04 


1.85- 


-2.20 


11.78 


2.72 




AzTEC.GSll 


1.69 


1.50- 


-1.95 


11.23 


4.37 




AzTEC.GS12 


4.50 


4.38- 


-4.70 


11.18 


3.08 


4.760 a 


AzTEC.GS13 


2.38 


2.25- 


-2.48 


11.16 


4.76 




AzTEC.GS16 


1.80 


1.55- 


-2.20 


11.04 


1.61 




AzTEC.GS17 


2.40 


2.20- 


-3.00 


11.38 


2.64 




AzTEC.GS18 


2.80 


2.45- 


-2.95 


11.09 


7.45 




AzTEC.GS19-l 


1.62 


1.27- 


-1.97 


11.18 


2.73 
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1.83 
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-2.10 
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0.63 
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LESSJ033219 
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-2.10 
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3.26 




LESSJ033243 


2.08 


1.90- 


-2.30 


11.53 


1.80 


2.122 f 



ture requires consideration of the average values. For complete- 
ness, in objects which required multiple-component fits we deter- 
mined separate average scalelengths and Sersic indices using re- 
sults from both the dominant component, and a single component 
fit to the whole system. We find that both of these methods pro- 
duce virtually identical results (average agreement within 5%), in- 
dicating that (sub)millimetre galaxy hosts are reasonably extended 
(( r i/2) = 4.3 ± 1.8 kpc; median r 1 / 2 = 4.1 kpc) disk galaxies 
((n) = 1.1 ± 0.4; median n = 1.0). These values were calcu- 
lated excluding the relatively low-redshift bulge-dominated sources 
AzTEC.GS27 and AzTEC.GS38 as explained above. If the two 
bulge-dominated sources are included, the conclusion obviously re- 
mains virtually unchanged with no change in scalelength, and only 
a slight increase in average Sersic index ((ri/ 2 ) = 4.3 ± 1.7kpc, 
(n) = 1.4 ±0.9). 

These new measurements of scale-length and Sersic index for 
our (sub)millimetre galaxy sample are in excellent agreement with 
those derived from ground-based if -band data by Targett et al. 
(2011). A further exploration of of the extent to which HST and 
ground-based imaging yield comparable results is presented in sub- 
section 5.1.1 below. 
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Figure 4. The Sersic index distribution for the 22 primary (sub)millimetre 
galaxy counterparts as derived from our two-dimensional modelling. The 
grey shading from darkest to lightest indicates i) the secure (i.e. P < 0.1) 
(sub)millimetre galaxy identifications, ii) the unsecure (i.e. P > 0.1) iden- 
tifications and, for comparison iii) the Sersic index distribution displayed 
by a complete mass-selected sample of 83 comparably massive galaxies at 
2 < z < 3, as deduced by Brace et al. (2012) from directly-comparable 
CANDELS iiieo imaging of the UDS field. 



Figure 5. The distribution of axial ratios for the 22 primary (sub)millimetre 
galaxy counterparts as derived from our two-dimensional modelling. 
The grey shading from darkest to lightest indicates i) the 17 secure 
(sub)millimetre galaxy identifications, ii) the 5 less secure (sub)millimetre 
galaxy identifications and, iii) for comparison, the axial-ratio distribution 
displayed by the 83 comparably-massive galaxies in the CANDELS UDS 
field at 2 < z < 3, as deduced by Brace et al. (2012), also using GALFIT 
applied to the directly-comparable CANDELS UDS Hieo imaging. 



The axial ratio distribution of our (sub)millimetre galaxies 
is plotted in Fig. 5. This provides an alternative (and model- 
independent) test of galaxy morphology. Although a detailed sta- 
tistical comparison is limited by small sample size, it can be seen 
that the distribution of axial ratios for the (sub)millimetre galaxies 
is very broad and flat, extending to values ~ 0.2, as expected for 
a population of disk galaxies (Sandage et al. 1970). In contrast, the 
axial ratio distribution for a mass-selected sample of galaxies with 
M* > 10 11 M Q in the redshift range 2 < z < 3 from Bruce et 
al. (2012) also plotted in Fig. 2, is peaked around (b/a) pea k — 0.7 
with tails extending to (b/a) — 0.1 and (b/a) = 1.0. 

Thus, while there has been much debate in the literature con- 
cerning why the axial-ratio distribution of disk galaxies at z ~ 2 
appears more peaked than expected from a population of randomly- 
oriented thin disks (see, e.g., Bruce et al. 2012; Law et al. 2012a,b) 
the picture presented by the (sub)millimetre galaxies seems to be 
relatively straightforward; judged either by typical Sersic index or 
axial ratio distribution, (sub)millimetre galaxies are massive, star- 
forming disks. We stress that this conclusion stands irrespective of 
what the mechanism for triggering or fuelling the star-formation 
activity in these galaxies might be (a point discussed further be- 
low). 

5.1.1 Comparison to ground-based observations 

In conjunction with the CANDELS HST imaging, deep ground- 
based A's-band imaging of GOODS-South has recently been ob- 
tained with the HAWK-I imager on the VLT. As described in sec- 
tion 2.3.3, these data were taken in excellent seeing (FWHM ~ 



0.45 arcsec) and reach a detection limit of ~26.0 (AB mag, 5-cr, in 
a 1 arcsec diameter aperture), a depth comparable to that achieved 
by the CANDELS HST ff-band imaging. The HAWK-I data have 
been reduced and processed as described in Fontana et al. (2012). 

This new deep HAWK-I imaging allows us to make a direct 
test of the extent to which the basic morphological parameters de- 
rived from modelling the high-resolution WFC3/IR i3i6o imaging 
agree or disagree from those that would be derived with GALFIT 
from poorer-resolution ground-based imaging available (albeit at 
slightly longer wavelengths). We have thus fitted axi-symmetric 
models to the K s -band images of the fourteen (sub)millimetre 
galaxies in our sample which lie within the HAWK-I imaging. The 
lower resolution of the ground-based data does not allow multiple 
component fits to the more complex systems, and therefore might 
be expected to be more comparable to the single-component model 
fits (see Table 2) attempted to the CANDELS data. 

The results of this 7\" s -band modelling are summarised in Ta- 
ble 4. The modelling was carried out in the same manner as de- 
scribed for the modelling of the ffieo imaging in Section 4. We 
find that the average half-light radii and Sersic index of both the 
single and dominant-component iJi 6 o-band model fits match those 
derived from the K s -band images to within 5%. A comparison of 
the HST and ground-based results for individual objects is shown 
in Fig. 6. These results suggest that, despite the extra detail re- 
vealed by the HST imaging, the basic sizes and morphologies of 
(sub)millimetre galaxies can be successfully recovered from the 
very best i\~ s -band imaging, explaining the similarity between the 
results presented here and those reported by Targett et al. (2011), 
and consistent with our finding that single-component fits to the 



© 0000 RAS, MNRAS 000, 000-000 



12 T.A. Targett, et al. 



Table 4. Results from two-dimensional modelling of AzTEC and LABOCA 
selected (sub)millimetre galaxies, as derived from the VLT HAWK-I K s - 
band imaging. Column 1 gives the source name. Column 2 lists the total 
host-integrated ii" s -band mag (AB). Column 3 lists the semi-major axis 
scalelength (half-light radius) of the host galaxy fit in kiloparsec. Column 4 
lists the best-fitting value of the Sersic index (n). Column 5 gives the axial 
ratio of the host galaxy. Column 6 gives the value of reduced xi for each 
model fit. 



Name K-mag r 1 / 2 Sersic n Axial ratio xi 

0.914 
1.036 
0.999 
0.787 
0.787 
1.298 
0.967 
1.068 
1.021 
1.082 
1.077 
0.858 
0.734 
6.542 



high-resolution CANDELS images correctly recovers the basic 
morphology of the mass-dominant disk-like component. 



5.1.2 Comparison to ULIRGs in CANDELS 

Recently, Kartaltepe et al. (2012) have studied the morphologies 
of a sample of 52 Ultra Luminous Infrared Galaxies (ULIRGs) in 
the CANDELS GOODS-South field, selected at 100 and 160/um 
from Herschel PACS imaging. Based primary on visual classifi- 
cation, Kartaltepe et al. (2012) reported that the majority of the 
GOODS-South ULIRGs were mergers or interacting sources, with 
a relatively small fraction of isolated disks. 

However, it is not clear that these results are in conflict with 
those found here for the (sub)millimetre-selected galaxies, and in 
fact it can be argued they are perfectly consistent. First, there is rel- 
atively little overlap between the two samples; the ULIRGs were 
selected at much shorter wavelengths and the resulting sample is 
therefore inevitably dominated by hotter and/or lower-redshift ob- 
jects, and such a sample may well contain a higher fraction of com- 
pact merger-driven starbursts and/or more significant bulges. There 
are in fact only six objects in common (AzTEC.GS06, GS08, GS 17, 
GS21, GS24 and LESSJ033219), and Kartaltepe et al. (2012) clas- 
sify three of these as isolated disks, and three as irregular, a re- 
sult very comparable to our own findings (albeit based on only 
a small subsample of our sources). Second, it is clear that visual 
classification inevitably places a stronger emphasis on high sur- 
face brightness features and morphological peculiarities than the 
multi-component modelling approach adopted here; several of our 
massive disk galaxies (especially those with clumps and/or com- 
panions) might well be classified as irregular and/or disturbed via 
purely visual classification. Third, Kartaltepe et al. (2012) did also 
use GALFIT to determine scalelengths and Sersic indices for their 
ULIRGs, and their results are not really very different from our 
own; they report a median R e — 3.3 kpc, making the ULIRGs only 



~ 20% more compact than our (sub)millimetre galaxies (median 
R e = 4.0 kpc), and a median Sersic index n = 1.4 (not as low as 
our own median Sersic n — 1.1, but still consistent with a disk- 
dominated population). 

5.2 Redshifts and stellar masses 

The redshift distribution of our (sub)millimetre galaxy sample 
is plotted in Fig. 7. We emphasize that ~ 50% of these red- 
shifts are spectroscopic, (including the one galaxy at z > 4.5, 
AzTEC. GS 12) and that Fig. 3 provides confidence that the pho- 
tometric redshifts are among the most accurate ever derived for 
(sub)millimetre galaxies. This redshift distribution is thus essen- 
tially complete and unbiased, and is broadly as expected given stud- 
ies of other, larger, (sub)millimetre-selected samples of galaxies 
(e.g. Chapman et al. 2005; Wardlow et al. 2011; Michalowksi et 
al. 2012b). Viewed together, Figs 4, 5, and 7 provide compelling 
evidence that (sub)millimetre selection primarily picks out star- 
forming disk galaxies at the redshift corresponding to the peak 
epoch of star-formation density in cosmic history. 

These disks are also all massive. As tabulated in Table 3, our 
SED fitting, assuming a Chabrier (2003) IMF, reveals all of the 
(sub)millimetre galaxies to have stellar masses M* > 10 11 M , 
with an average stellar mass (M*) = 2.2 x 10 11 ± 0.3 M . Again 
these results are consistent with much of the recent literature (e.g. 
Borys et al. 2005; Dye et al. 2008; Targett et al. 201 1; Michalowski 
et al. 2012a); the reader is referred to Michalowski et al. (2012a) 
for a detailed discussion of how the the choice of stellar population 
model, star-formation history and IMF can influence the estimation 
of stellar mass (the key points are that adoption of a Salpeter (1955) 
IMF increases stellar mass by a factor ~ 2, while adoption of a 
single-component star-formation history generally reduces stellar 
mass by a comparable factor). As discussed in the next sub-section, 
the stellar masses and sizes of the (sub)millimetre galaxies appear 
to be entirely plausible in the context of mass-selected galaxy sam- 
ples at comparable redshifts. 

5.3 The size-mass relation 

The average half-light radius and stellar mass of our 
(sub)millimetre galaxies is plotted in Fig. 9, overlaid on the 
distribution displayed by the general population of massive disk- 
dominated galaxies at 2 < z < 3 as derived by Bruce et al. (2012). 
The curves on this plot indicate the size-mass relations observed 
for early- and late-type SDSS galaxies at low-redshifts from Shen 
et al. (2003) (see figure caption for details). 

The Bruce et al. (2012) results illustrate the generally observed 
evolution towards smaller sizes at a given mass with increasing red- 
shift; massive disk galaxies at z ~ 2 seem to range in size from ob- 
jects which are still (just) consistent with the local relation to galax- 
ies almost an order of magnitude more compact. Viewed in this 
context, the (sub)millimetre galaxies appear reasonably extended 
for their mass ((ri/ 2 ) = 4.5 ± 2.5 kpc; median r 1 / 2 = 4.0 kpc), 
consistent with the sizes of the larger objects in the Bruce et al. 
(2012) sample, and indeed with other massive star-forming disks at 
z ~ 2 (Law et al. 2012a). In fact, the (sub)millimetre galaxies are, 
on average, only a factor < 1.5 smaller than present-day galaxies 
of comparable stellar mass. 

Clearly these results are at least consistent with the view that 
that most (sub)millimetre galaxies are simply drawn from the high- 
mass end of the normal star-forming disk galaxy population at 



AzTEC.GS16 


23.11 


5.1 


1.5 


0.62 


AzTEC.GS17 


21.89 


1.4 


1.7 


0.85 


AzTEC.GS18 


23.28 


3.1 


1.0 


0.70 


AzTEC.GS19-l 


21.85 


2.3 


1.5 


0.66 


AzTEC.GS19-2 


21.23 


3.7 


0.5 


0.65 


AzTEC.GS21 


21.37 


3.4 


1.4 


0.73 


AzTEC.GS22 


22.65 


3.7 


1.2 


0.40 


AzTEC.GS24 


24.07 


1.5 


0.8 


0.50 


AzTEC.GS26 


21.61 


8.2 


1.7 


0.36 


AzTEC.GS27 


20.05 


2.1 


3.5 


0.87 


AzTEC.GS30 


20.77 


4.5 


1.3 


0.70 


AzTEC.GS34 


21.11 


1.4 


1.2 


0.81 


AzTEC.GS35 


21.94 


3.3 


0.7 


0.74 


AzTEC.GS38 


18.32 


4.1 


4.6 


0.99 



© 0000 RAS, MNRAS 000, 000-000 



The properties of (sub)millimetre galaxies revealed by CANDELS 13 




2 3 
Sersic index (H-bond) / n 




2 4 6 8 

ho If — light rodii (H-bond) / kiloporsec 



Figure 6. Comparison of two-dimensional model Sersic index (left) and semi-major half-light radii (right) from (sub)millimetre galaxies in the Tiigo-band 
HST CANDELS data and the i^ s -band HAWK-I VLT imaging. The values shown for the WFC3/IR imaging represent those of the single-component fits to 
the (sub)millimetre galaxies, although the results remain unchanged if the dominant model component is adopted. 



z ~ 2 (as proposed, theoretically, by Dave et al. 2010). We now 
briefly consider whether our imaging provides any significant evi- 
dence for or against the viewpoint that the high star-formation rates 
displayed by the (sub)millimetre galaxies are the result of (or are 
significantly enhanced by) major galaxy-galaxy interactions (as, for 
example, inferred from the theoretical work of Narayanan et al. 
2010a,b, and arguably supported by the kinematical observational 
studies of Tacconi et al. 2008; Engel et al. 2010; Alaghband-Zadeh 
et al. 2012). 



5.4 Clumps, companions, and interactions? 

Our results show, beyond reasonable doubt, that (sub)millimetre 
galaxies are large, massive star-forming disk galaxies at z ~ 1 — 3. 
However, at least in principle, this does not necessarily prove that 
they all lie on the so-called "main-sequence" of star-forming galax- 
ies at z ~ 2 (e.g. Daddi et al. 2007; Elbaz et al. 201 1; Wuyts et al. 
2011). Even if essentially all (sub)millimetre galaxies are massive 
disks, at least some of their star-formation activity could be ex- 
ternally triggered. Indeed, as argued by Hayward et al. (2012) the 
(sub)millimetre galaxy population could contain a mix of objects 
including galaxies at the high-mass end of the star-forming main- 
sequence, along with other objects whose star-formation may have 
received a short-term boost from a galaxy-galaxy interaction (an, 
in addition, some straightforward source blends). 

The issue of galaxy interactions is arguably best addressed 
with the velocity information provided by spectroscopy (although 
even then the interpretation of the data at this epoch of generally 
violent star-formation activity is not straightforward). However, 
we can at least attempt to quantify the fraction of potential major 
galaxy mergers in our sample from our multi-wavelength broad- 
band imaging data, and our multi-component Sersic modelling. 

One simple way to do this is to investigate the relative ff-band 




Figure 7. The redshift distribution for our full sample of (sub)millimetre 
galaxies. The secure (i.e. P < 0.1) identifications are shown with dark 
grey shading, while the less secure (i.e. P > 0.1) identifications are indi- 
cated by light grey shading. 10 of these redshifts are the result of optical 
spectroscopy as summarized in Table 3, with the remainder being based on 
the SED modelling of the optical-near-infrared photometry as described in 
Section 4.2 (the reliability of which is supported by Fig. 3). 



flux density ratio of the primary and brightest secondary compo- 
nents seen in the WFC3/IR imaging of each object. This is shown 
(plotted in terms of magnitude difference) in Fig. 9. Here it can 
be seen that only ~ 1/3 of the (sub)millimetre disk galaxies have 
even a potential companion object which is more than ~ 0.1 times 
as bright as the dominant component. Treating if -band brightness 
as the best proxy we have for the stellar mass of individual sub- 
components, and bearing in mind that there is no evidence that 
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Figure 8. The location of our (sub)millimetre galaxy sample on the size- 
mass plane is indicated by the red star (with the error bar indicating the 1-cr 
range of the full size distribution). The black circles indicate the positions 
of the massive disk galaxies at z > 2 studied by Bruce et al. (2012). The 
solid blue and red lines show the size-mass relations deduced by Shen et al. 
(2003) for local early- and late-type galaxies respectively, with the dashed 
blue lines representing the 1-<t scatter for the late-type galaxies. 



all the brighter companion objects actually lie at the same red- 
shifts as the dominant (sub)millimetre disk-galaxy counterparts, 
it seems clear that only a minority of the (sub)millimetre galax- 
ies (probably < 25%) can be viewed as potentially involved in a 
major (i.e. > 1 : 3 mass ratio) galaxy-galaxy merger. This pro- 
vides further support for the conclusion reached by Michalowski et 
al. (2012a), based on the independent calculation of specific star- 
formation rates, that most (sub)millimetre galaxies are not primar- 
ily the result of extreme merger-driven starbursts, but in fact lie 
on the main-sequence of star-forming galaxies, with star-formation 
rates broadly as expected at z ~ 2 — 3 given their redshifts and 
stellar masses. 

Nevertheless, it is clear that, in HST imaging of this depth and 
quality, secondary components can be uncovered at some level for 
most of the galaxies in our sample. The key question is whether 
they are generally significant, or indeed unusual at this epoch of 
maximum star-formation activity in the gas-rich young universe. 

We can attempt to investigate this issue in two ways, first by 
reviewing the results of our attempts to model any apparently sig- 
nificant multiple-component structures in the ffieo imaging (as il- 
lustrated in Fig. 1, and shown for all objects in the Appendix in 
Fig. Al), and second by considering the multi-frequency imaging 
of each object (presented in Fig. A2). 

As discussed in Section 4.1, we attempted to model all appar- 
ently significant substructures or companions using variable Sersic- 
index galaxy models. This was done both to help obtain satisfac- 
tory model fits to the primary galaxy identification, and to see if 
any useful information could be gleaned from the modelling of the 
secondary sources. As can be seen from the summary of the mod- 



eling presented in Table 2, in practice models were fitted to sec- 
ondary components for approximately half of the (sub)millimetre 
sources. It can also be seen that the modelling was largely con- 
fined to secondary sources/components with Hu,o < 24.5. This is 
simply because the modelling of fainter subcomponents was nei- 
ther necessary nor practical (as a consequence of this it transpires 
that all 7 of the brighter secondary sources in the left-hand peak 
of Fig. 9 were modelled, while only 2 of the relatively fainter sec- 
ondary sources were modelled, with the flux-densities of most of 
the fainter clumps/companions being derived from aperture pho- 
tometry). An interesting feature of the results presented in Table 
2 is that the modelling of many of the fainter secondary sources 
delivered very low values of the Sersic index n < 0.5. A good 
example is the case of LESSJ033243, already highlighted in the 
second row of plots in Fig. 1 . Here the brightest secondary source 
appears to be a disk galaxy with n — 0.8, while the two even 
fainter objects (which to be distorted clumps), have n = 0.4. As 
a further example, the two subsidiary components of AzTEC.GS21 
have n = 0.3 and n = 0.2. By contrast, no primary galaxy iden- 
tification in the sample has n < 0.8. This suggests that Sersic in- 
dex can possibly be used to differentiate between genuine galaxy 
companions, and distorted clumps of excess emission. Some of the 
latter may be dwarf irregular galaxies, but they could equally just 
be viewed as "weather". If we apply the criterion that any gen- 
uine candidate for a major galaxy-galaxy interaction (or indeed sig- 
nificant (sub)millimetre source blending) must have a secondary 
component with a "sensible" Sersic index, then we are left with 
only 6 examples of disk galaxies with potential galaxy compan- 
ions (AzTEC.GS13, AzTEC.GS19, AzTEC.GS23, AzTEC.GS25, 
LESSJ033219, LESSJ033243), again representing ~ 25% of the 
sample. 

Second, moving away from the detailed modelling, a review 
of the multi-wavelength imaging presented in Fig. A2 reinforces 
the view that most complex substructure is blue, becoming less 
and less evident at longer wavelengths as the mass-dominant pri- 
mary galaxy identification becomes progressively more important. 
Thus, while at rest-frame UV wavelengths many of the galaxies in 
our sample look complex, clumpy and sometimes distorted, it is 
clear that much of this structure must involve relatively insignif- 
icant fractions of the stellar mass of the galaxy which hosts the 
(sub)millimetre emission. This sort of behaviour is well illustrated 
by AzTEC.GS08 (the first example shown in Fig. 1) where the 
multi-wavelength imaging shown in Fig. A2 shows that the clumps 
in the galaxy (which are too faint in Hieo to merit individual mod- 
elling) are very blue; viewed in the optical AzTEC.GS08 would 
certainly not be reported as a single isolated disk galaxy, and even 
at J125 it could easily be mis-interpreted as a merger or disturbed 
system. Indeed, if the signahnoise in our CANDELS WFC3/IR H- 
band imaging is artificially degraded to the level typically achieved 
with NICMOS (e.g. Swinbank et al. 2010) it becomes clear that, 
even at if -band, imaging of inadequate depth would fail to prop- 
erly uncover the underlying mass-dominant disk galaxy, instead 
just revealing the highest surface-brightness features. This is per- 
haps not surprising; at z ~ 2.5 the H-band is sampling the rest- 
frame B-band, which is why it would really be better to move to 
K or ideally IRAC wavelengths to determine the morphology of 
the mass-dominant galaxy, if only imaging of the necessary angu- 
lar resolution were available at these wavelengths (deep ~ 4 /im 
imaging with the James Webb Space Telescope will therefore be 
ideal for addressing this issue). 

We close by noting that the prevalence of clumpy/complex 
structures in star-forming galaxies at z ~ 2 is certainly not con- 
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Magnitude difference (H-band) 

Figure 9. The magnitude difference between the brightest and second- 
brightest //-band components for all 20 disk-dominated (sub)millimetre 
galaxies. As in previous figures, for completeness, the lighter grey shading 
is used to indicate the less reliable galaxy identifications (although clearly 
the basic result shown here is unaffected by this minor uncertainty). Only 
~ 1/3 of the (sub)millimetre galaxies have even a potential companion 
galaxy contributing > 10% of the //-band flux density of the primary mass- 
dominant disk galaxy. 

fined to (sub)millimetre galaxies. Indeed, as shown for example by 
Wuyts et al. (2012), most lower-mass star-forming galaxies at these 
epochs display clumpy, blue substructure, the origin and impor- 
tance of which remains the subject of vigorous theoretical debate 
(e.g. Ceverino et al. 2012; Genel et al. 2012). In this sense, while 
it can be argued that ~ 25% of the (sub)millimetre galaxies in our 
sample could be involved in a significant galaxy-galaxy interaction, 
in general the galaxies selected at (sub)millimetre wavelengths do 
not seem unusual for star-forming galaxies at z ~ 2 — 3 (except 
that they are among the most massive and largest galaxies at that 
epoch). There is in fact now considerable doubt over whether ma- 
jor mergers are common enough at z ~ 2 — 3 to explain the preva- 
lence of (sub)millimetre galaxies at this epoch (e.g. Hopkins et al. 
2010). Moreover, as pointed out by Cen et al. (2012), that even for 
those z ~ 2 (sub)millimetre galaxies which arguably do show some 
evidence of companions or tidal interactions, it does not necessar- 
ily follow that interactions and/or mergers are the primary cause 
of their high-starformation rates. Against the backdrop of our rela- 
tively peaceful present-day Universe, it is clear that galaxy-galaxy 
mergers can have a dramatic effect (e.g. in the production of the 
low-redshift ULIRGs first uncovered by the IRAS satellite), but 
amid the general frenzy of star-forming activity in the gas-rich uni- 
verse at z ~ 2 (e.g. Dekel et al. 2009) the relative impact of merg- 
ers on the star-formation rate of massive galaxies may be largely 
inconsequential. 



Spitzer (IRAC 8 fim) imaging to identify galaxy counterparts for 
the (sub)millimetre galaxies, we have then utilised the new CAN- 
DELS WFC3/IR imaging in two ways. 

First, the addition of new deep near-infrared photometry from 
both HST and (at K-band) HAWK-I on the VLT to the exist- 
ing extensive GOODS-South multi-frequency database has enabled 
us to derive improved photometric redshifts and stellar masses 
for all the (sub)millimetre sources. Our results confirm that the 
(sub)millimetre sources are massive ((M*> = 2.2x 10 11 ±0.3 M ) 
galaxies at z ~ 1 — 3. 

Second, we have exploited the depth and superior angular 
resolution of the WFC3/IR i/ieo imaging to determine the sizes 
and morphologies of the galaxies at rest-frame optical wavelengths 
Arcst > 4000A. Specifically, we have fitted two-dimensional axi- 
symmetric galaxy models to the WFC3/IR images, varying lu- 
minosity, axial ratio, half-light radius ri/2, and Sersic index n. 
Crucially, the wavelength and depth of the WFC3/IR imaging en- 
ables modelling of the mass-dominant galaxy, rather than the blue 
high-surface brightness features which often dominate optical (rest- 
frame ultraviolet) images of sub-mm galaxies, and can confuse vi- 
sual morphological classification. As a result of this analysis we 
find that > 95% of the rest-frame optical light in almost all of the 
(sub)millimetre galaxies is well-described by either a single expo- 
nential disk (n ~ 1), or a multiple-component system in which the 
dominant constituent is disk-like. We demonstrate that this conclu- 
sion is completely consistent with the results of recent high-quality 
ground-based _R"-band imaging studies of (sub)millimetre galaxies 
sampling even longer rest-frame wavelengths (Targett et al. 201 1), 
and explain why. We also briefly compare our findings with the 
results of recent morphological studies of high-redshift ULIRGs 
selected at much shorter (~ 10 x) far-infrared wavelengths (Kar- 
taltepe et al. 2012). 

The massive disk galaxies which host strong (sub)millimetre 
emission are reasonably extended ((ri/2) = 4.5 ± 2.5 kpc; me- 
dian ri/2 = 4.0 kpc), consistent with the sizes of other massive 
star-forming disks at z ~ 2. In many cases we find evidence of 
blue clumps within the sources, with the mass-dominant disk com- 
ponent becoming more significant at longer wavelengths. Interest- 
ingly, the only two bulge-dominated sub-mm galaxies are also the 
two lowest-redshift sources in the sample (z ~ 1), a result which 
may reflect the structural evolution of high-mass galaxies in gen- 
eral. 

Our main result that (sub)millimetre galaxies are almost uni- 
versally massive, star-forming disk galaxies does not exclude the 
possibility that some fraction of these objects have their star- 
formation boosted by interactions/mergers. However, our imaging 
data indicate that only a minority (< 25%) of the (sub)millimetre 
sources in the GOODS-South sample could be regarded as involved 
in a major galaxy-galaxy interaction. Taken together, our results 
support the view that most sub-mm galaxies at z ~ 2 are simply 
the most extreme examples of normal star-forming galaxies at that 
era. 



6 CONCLUSIONS 

We have exploited the Hubble Space Telescope {HST) CAN- 
DELS J and f/-band WFC3/IR imaging to study the properties 
of (sub)millimetre galaxies in the GOODS-South field which have 
been revealed by the AzTEC 1.1 mm and LABOCA 870 /im sur- 
veys of the Extanded Chandra Deep Field South (Scott et al. 2010; 
Weiss et al. 2009). After using the deep radio (VLA 1.4 GHz) and 
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APPENDIX A: IMAGES AND GALAXY MODEL FITS 

In this Appendix we first show the CANDELS Hieo images of all 
the galaxies in our (sub)millimetre sample, along with the model 
fits and residual images (Fig. Al), and then also provide a multi- 
wavelength montage for each object to reveal the wavelength- 
dependence of the structure seen in the imaging (Fig. A2). 
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A -2 -A A -2 -A A -2 -A 

Figure Al. Two-dimensional modelling of the AzTEC and LABOCA selected (sub)millimetre galaxies in CANDELS GOODS-South. The left-hand panel 
shows the WFC3/IR //igo-band image centred on the (sub)millimetre galaxy counterpart. The centre panel shows the best-fitting two-dimensional galaxy 
model (convolved with the HSTVSY. The right-hand panel shows the residual image after subtraction of the model from the data. All panels are 12.0" X 12.0", 
and the images are shown with a linear greyscale in which black corresponds to 10-ct above, and white to l-a below the median sky background value. 
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4 -2 -A A -2 -A A -2 -4 

Figure Al. - continued 
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Figure Al. - continued 
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